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ABSTRACT 


In this project, the study has been carried out to investigate the effects of turbine 
runner parameters on the design. For determining the parameters effects a Francis 
turbine runner for high specific speed has been designed and the parameter has been 
studied and numerically analyzed. Design of turbine runner blade is done using Bovet 
approach. Blade geometry of runner is generated using ANSYS BladeGen software. 
Numerical analysis has been carried out to determine different parameter's effects on 
turbine performance. By varying blade outlet angle efficiency and sediment erosion 
has been compared to design an efficient model of Francis turbine — runner for high 


specific speed. 
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СНАРТЕК 1 INTRODUCTION 


1.1 Context 


Energy demand of the world is increasing day by day and the use of renewable energy 
sources as hydropower can be a option for the fulfillment of energy demands. As 
hydroelectricity is the generation of electric power by use of energy stored in water. 
The main components of a hydropower is the turbine. Generally the basic types of 
turbines are : pelton, francis, crossflow, kaplan , turgo,etc. The selection of the use of 
turbines is generally site dependent and according to the turbine selection chart francis 


turbine has wide range of selection criteria. 
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Figure 1 Turbine Selection Chart [2] 


The variation range of selection of Francis turbine according to head, flow rate and 
power output is a better option for easy and efficient type of hydropower construction. 
A Francis turbine consists of mainly following components : 

Spiral casing: The spiral casing around the runner of the turbine 1s known as the volute 
casing or scroll case. Throughout its length, it has numerous openings at regular 
intervals to allow the working fluid to impinge on the blades of the runner. These 
openings convert the pressure energy of the fluid into momentum energy just before 
the fluid impinges on the blades. This maintains a constant flow rate despite the fact 
that numerous openings have been provided for the fluid to enter the blades, as the 
cross-sectional area of this casing decreases uniformly along the circumference. 

Guide or stay vanes: The primary function of the guide or stay vanes 1s to convert the 
pressure energy of the fluid into the momentum energy. It also serves to direct the 
flow at design angles to the runner blades. 

Runner blades: Runner blades are the heart of any turbine. These are the centers where 
the fluid strikes and the tangential force of the impact causes the shaft of the turbine to 
rotate, producing torque. Close attention in design of blade angles at inlet and outlet is 
necessary, as these are major parameters affecting power production. 

Draft tube: The draft tube is a conduit that connects the runner exit to the tail race 
where the water is discharged from the turbine. Its primary function is to reduce the 
velocity of discharged water to minimize the loss of kinetic energy at the outlet. This 
permits the turbine to be set above the tail water without appreciable drop of available 


head. 


1.2 Background 


since the world is giving priority to renewable energy sources rather than non- 
renewable fossil fuels. Hydropower being the most important energy source to meet 
the needs and criteria of energy in the present world .Francis turbine in comparison to 
others , due to performance, operating expenses and efficiency range allows more 
versatility in choice of head, flow rate and power. In context of Nepal most of the 


francis turbines used are for low specific speed and researches related to high specific 


speed 15 not available. This deficiency in efficient design capabilities has led to more 
research and project to expand design techniques for high specific speed francis 
turbine. 

Performance and efficiency wise it has been seen that francis turbine for low head has 
more advantages but due to the lack of proper researches and design techniques of 
high specific speed francis runner it's use has been limited so this project entitled 
"DESIGN OF HIGH SPECIFIC SPEED FRANCIS TURBINE FOR LOW HEAD 
APPLICATIONS"' can develop an efficient design of high specific speed francis 


runner to increase it's use ш micro hydro. 


1.3 Objectives 


1. Design low head Francis turbine. 
2 .Numerical analysis of the runner using CFD technique. 


3. Compare efficiency and sediment erosion rate by varying blade outlet angle. 


1.4 Significance/Scope 


1. Preliminary design will be done using Bovet approach. 
2. Profile generation will be done using BladeGen in ANSYS. 


3. Design of runner below 50m head. 


1.5 Limitations 


1. Among different components of Francis turbine runner will be designed. 

2. Parametric designs affecting runner design only will be considered. 

3. Runner design will be for high specific speed and using Bovet approach. 

4. There is no experimental validation of the results obtained from CFD. 

5. Economic analysis of cost of efficiency change with increase or decrease in erosion 


has not been done. 


CHAPTER 2 LITERATURE REVIEW 


2.1 Design approach 


In the design and development of Francis turbine various researches and works have 
been done to attain the desired efficient model of low specific speed Francis runner. 
some of the results and conclusions of the works carried out are as follows: 
Lars Froyd (2016) states important features in a low head design: 
1) The curved shape of the blade leading edge, which results in different inlet 
diameters at the top and the bottom of the turbine inlet. 
2) The small difference between inlet and outlet diameter giving a high shroud 
curvature. 
3) The straight vertical or sometimes even outward sloping shape of the shroud at 
the turbine outlet.[1 | 
Ayancik Fatma(2014) found parametric design based on numerical calculation 
indicates that hydraulic design parameters causes different effects on the blade shape 
and the turbine performance. Parametric relations between the beta and theta angle 
shows that the beta angle distribution has an important role in the design process. [2] 
Kazayuki Nakamura and Sadao Kurosawa (2009) optimized runner shape 
automatically enables on the constraints of the trade off relationship between the 
runner loss and the draft tube loss. The runner has great improvement on the turbine 
operating efficiency in the wide operating ranges.[3] 
From these literature review it is clear that yvarious works have been done in the field 
of francis turbine and runner being the most complex and affecting component needs 


more researches and works based on the high specific speed runner is yet to be done. 


MERIDIONAL PLANE 
In the first step of classical design method of a Francis type hydraulic turbine, by the 
assumption of axial symmetry, the meridional view is obtained. The important 


dimensions of the runner are calculated by the BOVET. After the dimensions are 


determined, we can find the entrance and exit limit curves, that in our case the 
BOVET method has been used. 

Fundamental parameters in turbine design are: 

H: Hydraulic head [m] 

Q: Turbine flow rate [m3/sec] 

N: Turbine turns [rpm] 


Regarding these parameters, the dimensionless specific speed n, can be defined as 





follows: 
| N(P)1/2 
ET е 


Specific rotational speed, ‘n’? is used to calculate the values of runner channel 
parameters in the meridional plane. Bovet approach is capable to calculate runner 
channel parameters in dimensionless form as function of ‘n’. Method uses specific 
radius of 'Roe'which is the distance between rotational axis of the runner and 
intersection of shroud and trailing edge of the blade. Specific radius can be calculated 


from equation (2). 


R= (2 = ) (2) 


п Маге 





Where, voe represents specific volumetric flow rate. It is taken as 0.270 for this study. 
To obtain dimensional results, parameters of meridional channel are multiplied with 
specific radius of the runner. By this way whole variables of the runner channel are 
reduced to functions of net head, flow rate and rotational speed. 


Table 1 Equation of meridional plane parameters [3] 


bo = 0.8(2-n)n Toi7ymi70.7 + 0.16/(n+0.08) 





li=3.2 + 3.2(2-n)n 1-2.4-1.9(2-п)п 


The correlation between specific rotational speed and meridional channel parameters 


are not linear as seen in Table 1. So this kind of relation between 'n' and Ше 
parameters causes logarithmic variations and naturally this leads to possibility of 


designing different types of blades. The BOVET well known relations give the 


geometry of hub and shroud curves in Ше runner blade zone. The form of limit curves 


is defined by equation: 


2--308(1-2) (1-5) (4) 
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Figure 2 Important dimensions of runner channel [2] 


After these necessary calculations are completed, runner channel meridional plane can 
be created. Then, intersection points of runner blade leading and trailing edges with 
hub and shroud were determined, respectively. Therefore, runner blade meridional 
profile was placed on runner channel meridional plane. It is known that the most 
important design parameter of a runner blade is profile throughout leading to trailing 
edges. For this purpose, Bovet separates the runner blade with streamlines which have 
equal flow rate between each other and this is main characteristic of these streamlines. 
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Seven streamlines were used which are named from hub to shroud respectively ii, Т”, 


66299 
l 


m, e’’, e’, e where “i” spline represents hub curve while “e” curve representing shroud 
(Fig. 4). It guides to design more sensitive geometry while constructing runner in 
CAD software. Completing construction of meridional profile, it is also need a 
determination for designate the angle of specific points on the blade. Recent studies 
are showed that the angle of leading and trailing edges are highly significant when 
designing high efficient runner. The design approach calculates radiuses and angles at 


intersection points of streamlines with trailing and leading edges by using velocity 


triangles. 
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Figure 4 Streamlines on channels' meridional plan [2] 


2.2 Francis turbine terminology 
There are number of terminology related to the Francis turbine. 
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Figure 5 Bladegen generated meridional profiles 


Hub: It is the side of axis of turbine runner where shaft is attached. 

Shroud: It is the opposite of the hub as shown in above figure. 

Leading edge: It is the edge of the Francis runner blade at inlet through which water 
strike the blade. 

Trailing edge: It is edge of the Francis runner blade through water discharges to the 
draft tube. 

Pressure side: It is the side of the blade at which incoming water strike through inlet. 
Suction side: It is the side opposite to the pressure side. 

Meridional view: The view of the blade in the radial and axial coordinate system. 
Blade span: It is the length of the blade at the inlet through which water impinges at 
blade. 


2.3 Computational Fluid Dynamics 
Computational Fluid Dynamics (CFD) is one of Ше branches of fluid mechanics that 


uses numerical methods and algorithms to solve and analyze problems that involve 
fluid flows. Millions of calculations required to simulate the interaction of fluids and 
gases with the complex surfaces studied in science and engineering by iterative 
process. However, even with simplified equations and high speed supercomputers, 


only approximate solutions can be achieved in many cases. (Khanna, 2011) 


2.3.1 CFD Structure 
СЕО 15 structured around the numerical algorithms used to tackle fluid problems. The 
CFD software includes sophisticated Graphics User Interfaces (GUI) to input problem 
parameters and to examine the results. The basic components are: 

е  Pre-processing 

е Solver 


е Post processing 


2.3.1.1 Pre-processing 
It 1s the first step of analyzing the problem. The activities performed at the pre- 
processing stage are: 
е Definition of the geometry of the region: the computational domain 
е Grid generation: the subdivision of the domain into a number of smaller, non- 
overlapping sub domains. 
e Definition of fluid properties 
e Specification of appropriate boundary condition at cells, which may coincide 
with or touch the boundary 
2.3.1.2 Solver 
The main part of the CFD structure 1s the solver. FLUENT, CFX and POLYFLOW 
are some of the types of solvers. CFX solver 15 developed independently by ANSYS. 
It is control-volume based for high accuracy and rely heavily on a pressure based 


solution technique for broad applicability. The СЕХ solver uses finite elements (cell 


vertex numeric), similar to those used in mechanical analysis, to discretize the domain. 
The numerical solution of Navier-Stokes (NS) equations in СЕО usually implies а 
discretization method: it means that derivaties in partial differential equations are 
approximated by algebraic expressions which can be alternatively obtained by means 
of finite. difference of the finite element method. The result 1s a set of algebraic 
equations through which mass, momentum, and energy transport are predicted at 
discrete points in the domain. The governing equations are non-linear and coupled, 
several iterations of the solution loop must be performed before a converged solution 


15 obtained. 


2.3.1.3 Post- processing 
It is the final step in CFD analysis. The result obtained from the solver is interpreted to 


produce CFD images and animations. (Khanna, 2011) 


2.3.2 Discretization Methods in CFD 

There are three discretization methods in CFD: 
e Finite difference method (FDM) 
е Finite volume method (FVM) 
e Finite element method (FEM) 


2.3.2.1 Finite element method (FEM) 

А finite element method (FEM) discretization is based upon а piecewise 
representation of the solution in terms of specified basis functions. The discretization 
involves division of computational domain into smaller domains (finite elements) and 
the solution in each element is constructed from the basic functions. The actual 
equations that are solved typically obtained by restating the conservation equation in 
weak form: the field variables are written in terms of the basic functions; the equation 
is multiplied by appropriate test functions, and then integrated over an element. Again, 
a system of equations is obtained (usually for nodal values) that must be solved to 


obtain a solution. 
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2.3.3 Turbulence Models 

The CFD software include the package to model the fluid flow phenomena under the 
turbulent models. Usually turbulent numerical simulation consists of two main parts, 
namely: Direct Numerical Simulation (DNS) and Indirect Numerical Simulation 
(INS). DNS has a precise calculated result, but the whole range of spatial and temporal 
scales of the turbulence must be resolved which requires a very small time step size. 
So, this is not suitable for CFD simulations. There are three different types of 
simulated methods under the indirect Numerical Simulation which are large eddy 
simulation (LES), Reynolds-averaged Navier Stokes (RANS) and detached eddy 
simulation (DES). In order to simulate turbulent flows, theoretically, the 
computational domain should be big enough to contain the biggest eddy. Meanwhile, 
the mesh should be small enough to find out the smallest eddy. The equation of 
Reynolds-averaged Navier Stokes (RANS) is defined as: 


DU; OP , 9 | (ди 9UjV тт 
p Dt OK. Ox и (5: + 22) p ar; [3] 








It is the oldest and most common approach to turbulence modeling; the left hand side 
of the equation describes the change in mean momentum of fluid element and the right 


side of the equation is the assumption of mean body force and divergence stress. 


pul jul; is an unknown term and called Reynolds stresses. The RANS equation is not 
closed due to the presence of stress term, so it requires a turbulence model to produce 
a closed system of solvable equation. (ANSYS, 2009) 

The turbulence model selected 15 SST which is combine form of К-е and k-@ using Ше 
advantage of both model. 

2.3.3.1 Shear stress Transport (SST) k-@ Model 

The shear- stress transport (SST) k-@ model was developed by Menter to effectively 
blend the robust and accurate formulation of the k- @ model in the near —wall region 
with the free-stream independence of Ше К-в model in Ше far field. To achieve this, 
the k-e model is converted into a k-@ formulation. The SST k-@ model is similar to 


the standard k-@ model, but includes the following refinements: 
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е The standard k-@ model and the transformed К-е model are both multiplied by 
a blending function and both models are added together. The blending function 
is designed to be one in the near wall region, which activates the standard К-о 
model, and zero away from the surface, which activates the transformed К-е 
model. 

е The SST model incorporates a damped cross-diffusion derivative term in the o 
equation. 

е The definition of the turbulent viscosity is modified to account for the transport 
of the turbulent shear stress. 


e The modeling constants are different. 


2.3.4 Erosion model 


CFX solver has two inbuilt erosion model in it namely Finnie and Tabakoff model. 


2.3.4.1 Tabakoff and Grant Erosion Model 

In the erosion model of Tabakoff and Grant, the erosion rate E 1s determined from the 
following relations: 

E = К (у). М;“.сов (у) [1-R^r] + (Уру) 

Here, Е 15 the dimensionless mass (mass of егодед wall material divided by the mass 
of the particle), Vp is the particle impact velocity, y is the impact angle in radians 
between the approaching particle track and the wall, yo being the angle of maximum 
erosion, КІ to k4,k12 and yo are model constants and depend on the particle/well 


material combination. (ANSY 5,2006) 
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СНАРТЕК 3 METHODOLOGY 





Meriodional Plane design 
СЕО Analysis 


Expected 
outcome 














Final Design 


Figure 6 Flowchart of Francis Turbine Runner Blade Design 


The various Francis turbine set up at different hydropower stations for high head is 
being observed during the field visit and the runner design will be generated using 
Bovet approach. Based on the pre-existing designs and literature review design of the 
runner for high specific speed has been done in Auto CAD. The designed runner will 
be analyzed by varying the parametric constraints affecting the design. Based on the 
results obtained from numerical analysis an efficient design of Francis runner will be 


proposed. 
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Technigues used for efficient design of the runner are: 

Bovet approach: Bovet method uses empirical equations to obtain parameters of 
Francis type hydro turbine runner [4]. In Bovet method, the dimensionless specific 
speed value is the main parameter to determine whole turbine dimensions. 

BladeGen: BladeGen is a feature in ANSYS by which the single blade profile of the 
runner can be developed. Input parameters of head and flow rate are considered. Prior 
to the mesh generation for the design and CFD analyses, boundaries of the flow 
passage should be defined. BladeGen module of ANSYS was selected for the 
designing because of its comparative advantages. BladeGen provides rapid 
simulations, modification and optimization when used with ANSYS CFX. 

ANSYS BladeGen is a geometry creation tool that is specialized for turbo machinery 
blades. It is used to generate blades of radial turbine, axial, impeller, radial diffuser 


desired vane. 


Initial Meridional Configuration Dialog ? X 
Simple Axial Normal Axial Radial Impeller 
Radial Turbine Radial Diffuser Deswirl Vane 


Cancel 











А | 0.000000 


z- [180.000000 
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Figure 7 First Bladegen interface to input main dimensions of runner 
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Figure 8 Graph produced in BladeGen 


The high specific speed Francis runner has not been developed and the lack of 


researches has led to this project. 


3.1 Finding Main Dimensions 


Plots for Flow rate and head for different existing Francis turbine were created. And 
the region with maximum intensity was chosen to set the О and H of the design. 


Table 2 Design Parameters 


Specifications Value 


Volumetric flow rate (Q) 1.2 m/s 


Rotational speed 750 rpm 


Specific speed (n) 186 (from equation 1) 


Speed number 0.72 (Low head range) 


Specific volumetric flow rate (qoe) 


Specific radius (Кос) 0.31 (from equation 2) 
Meriodional ~— 0.6187 m 
channel parameters 
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0.964 m 


0.868 m 


These dimensions obtained from Bovet method has also been compared with Eulers 





method and after satisfactory result main dimensions are fixed. 


3.2 Modeling the Turbine 


The modeling of the turbine include the creating 2D view in the BladeGen software 
which is then converted into the 3D view by applying the beta distribution. This 
geometry is exported into the Turbogrid which is then meshed using Automotive 


Topology Meshing (ATM) method. 





Figure 9 3D geometry of runner 


3.3 Creating Domain and Mesh 


since very large number of mesh elements are needed to simulate the Francis turbine 
the domain with only one passage is created that has a blade in the middle. The 
domain starts from hub wall ends 1n shroud wall along the span wise direction while 
its starts from a small distance from the leading edge and ends in a small distance from 


the trailing edge. The calculations in the domain with whole turbine comprising wicket 
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gates, stay vanes, spiral casing and draft tube would increase Ше computational cost as 
it would take a huge amount of time and require a very powerful computer. 
Since the turbine is rotationally periodic it is assumed that the flow in the Francis 


turbine too is rotationally periodic and simulating in only one passage is enough to 


visualize the flow in the whole runner. 


ed 
2 





0 0.050 0.100 (m) 





0.025 0.075 


Figure 10 Mesh obtained from TurboGrid 


3.4 Problem Definition in CFD 
Problem definition refers to the setting of the analysis type, adding boundary 


conditions, choosing turbulence and erosion model in the CFD model. Following were 


chosen in ANSYS CFX to set up the problem. 
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3.4.1 Choosing Erosion Model 


In CFX, there are two namely: finnie and tabakoff model of erosion. The tabakoff 
model gives the better control on the variable of sediment and control structure. For 


our purpose the more control available will be better. So, tabakoff model is chosen. 
3.4.2 Choosing the Turbulence Model 


The standard k-e model is better for flow far from the boundary layer and the standard 
k-@ model is near the boundary layer. The Shear Stress Transport (SST) model utilize 
the positive aspect of the both the models. In Francis turbine runner the flow of water 
occurs both near the wall and far from the wall in the blade surface and the hub 


shroud surface. So, the shear stress transport (SST) model 1s chosen for the simulation. 
3.4.3. Boundary Conditions 


Analysis type: Steady State Analysis 


Domain: 

Fluid & particle definition Water, Quartz 
Reference pressure ] atm 

Erosion modeling Tabakoff erosion model 
Average diameter of quartz 100 micron 

Shape factor Off 

Turbulence model SST model 

Analysis type Steady state 

Drag force Schiller Naumann 
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Inlet boundary: 


Flow direction: Cylindrical components for sand 
No. of position 10000 
Mass flow rate of quartz 0.07 kg/sec 


Particle control: 


Number of integration steps 10 
Maximum particle Integration time step 10s 
Tracking time 108 
Tracking distance 10m 
Maximum number integration step 10000 





0 0.150 0.300 (m) 0 0.150 0.300 (m) 


0.075 0.225 0.075 0.225 





Figure 11 Boundary Conditions 


3.5 Running Simulation on CFX 
The above selected boundary conditions are used in the CFX pre, and .def 


file 1s created which 1s then fed into the solver. The CFX solver solve the 
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problem and to give the result Ше. The result thus, obtained is analyze ш 
ше СЕХ post to view Ше velocity distribution, pressure distribution, 


erosion, etc. as required. 


3.6 Observation of Trends of Erosion 


First the variation of erosion with change in blade outlet angle B2 was 
observed and a value of D» for minimum erosion was taken for the second 


kind of model creation. 
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CHAPTER 4 RESULT AND DISCUSSION 


4.1 Efficiency VS Blade Outlet Angle 


6 
Efficiency 





17.5 18 185 19 195 20 20.5 21 21.5 22 22.5 23 23.5 24 24.5 25 25.5 
Blade Angle В, 


Figure 12 Graph efficiency VS outlet angle 


By varying Ше outlet angle from 18° to 25° with 1^ span the efficiency of different 
modeled blade were obtained. The efficiency VS blade outlet angle was plotted and 
maximum efficiency of the blade was obtained at 23 blade outlet angle and maximum 
efficiency of 76%. 

The trend of efficiency variance with blade angle was: increase in blade angle resulted 
in the increase in efficiency from 18° to 23° which in the above graph is the line 
having positive slope. On further increasing the outlet blade angle from 23' the 
efficiency was found to decrease thus indicating maximum efficiency at blade outlet 


angle 23“. 
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4.2 Sediment erosion rate density at various blade outlet angles 
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Figure 13 Sediment erosion rate density comparison 


The above figure 11 shows the sediment erosion rate density on the suction side of the 
blades as the blade angle has been increased it has been found that the sediment 
erosion effect on the blade has decreased. Considering the efficiency and sediment 


erosion density rate optimum blade design at 23 outlet angle has been chosen. 
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CHAPTER 5 CONCLUSION AND RECOMMENDATION 


5.1 CONCLUSION 


With multiple simulations in different models and observing the values of erosion rate 
density in the blades and efficiency, relation has been figured out between p» with 
efficiency and erosion rate density. No clear relation of these parameters can be 
drawn. However, weak relations сап be drawn. Also clear and strong relation of Во 
with efficiency can be seen . 
The results can be concluded as follows: 
1) Bovet method is a method that uses empirical relations to determine 
meridional plane. The results obtained from CFD analysis are shown. 
2) Higher velocity in the suction sides cause sediment erosion in the suction sides 
of blades. 
3) Negative pressure zone in suction side of leading edge and also in trailing edge 
can cause cavitation. Cavitation formation must be considered 
4) The erosion was least with no significant decrease in efficiency in the blade of 


В2 value of 23°. 
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5.2 Recommendations 
1) Considerations of draft tube and spiral casing may give more accurate results and 


improved efficiency. 

2) Different erosion models can be considered to compare the sediment erosion 
density rate in the blades. 

3) Increasing the number of model blades can improve the understanding of the 


relations of D2 with efficiency and erosion rate density. 


24 


APPENDIX А 


3D MOELLING OF A FRANCIS RUNNER 


Designing of a Francis Runner with Hydrodynamic theory 


This theory is to determine the runner dimensions given the flow rate О and head H 


(Gogstad, 2012). 


Steps: 
е Choosing the outlet blade angle (between 13? to 24°) and peripheral velocity of 
blade inlet U» (between 35 to 42 m/s) 
139 <B2<22° lowest value for highest head 
35m/s <U2 <42 m/s Highest value for highest head 
е Flow at best efficiency point (Q*) = Q/1.2 [Q = 1.1 to 1.3 times Q*] 


For best efficiency point whirl component of outlet velocity is zero. С = 0 


U, =a, - R, 





Figure: Outlet velocity triangle 
е The meridional velocity at outlet Cm2 = U2 tan(B2) 


4Q* 


пСт2 





е The outlet diameter of the runner Р» = 


U2 





e The rpm of the runner = T 
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The rpm of the runner should be synchronized with the freguency of the national 


grid. For Nepal, it is 50 Hz. The frequency depends upon no. of pole pairs in Ше 


Number of pole pair (Z)*rpm 


generator and the rpm as f= я 


n(corrected) = 3000/2 
Accordingly the outlet diameter, outlet peripheral velocity, outlet meridional 
velocity have to be adjusted as: 


2400" 


D» (corrected) = m m 
z 





U» = (П D» n(corrected))/60 

См? (corrected) = U2 * tanpo 

The angular velocity of the runner 15 o= 2лп/60 

@*= o/ N2gH 

The speed number of the runner 15 Q = o NO 

The Net Positive Suction Head Required (NPSH) = 10 — hs where hs is the suction 
head or submergence height 


The empirical relation for the submergence height 1s: 


NPSH = 


(а.стг2+0.0,2) 
2 


if Q < 0.55 then a = 1.12, b=0.055 


If > 0.55 then a =1.12, b= 0.1.О“ 
The meridional velocity at the outlet is to Бе made 10% greater than that of the 


inlet. Hence we find Ше meridional velocity at Ше inlet as 





Ста“ Tm (10 % meridional acceleration) 


The reaction degree of Francis turbine 1s the ratio of pressure energy to the total 
energy available to the runner at the inlet. Usually the reaction degree lies in the 
periphery of 0.5 ie 50% K.E & 50 % pressure energy. The total energy per mass at 
the inlet is g * h and the kinetic energy per mass is Ci?/2. Hence the absolute 


velocity at the inlet 1s 


Cı =4(1 — reaction degree) x 2gH 


The velocity triangle at the outlet of the runner 1s shown: 
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Figure: Velocity triangle in Francis Turbine 


The whirl component of the velocity at Ше inlet is given by: 


Cul = V С — Сто“ 


е The Euler's equation is given as: 
Сасык HIED 
gH 
e Тһе whirl component at the outlet is zero Со = 0. Generally used value of 
efficiency is 0.96. The peripheral velocity of runner at inlet as: 
H 
e U=- 
Cu1 
e The diameter at the inlet 1s: 
_ 60U4 
Tm 


е 1 


e The height of the runner at Ше inlet 15: 
Q* 


TLD41Cm1 





e b= 
e The inlet angle of Ше blade is from the velocity triangle 


е В, = tan ! (е) 


Од Сил 
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Figure: Turbine main dimensions 
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Runner vanes 


Figure: Velocity triangles 
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APPENDIX В 


Table: Points plotted for limit curve 
о ff 
0.06 |0.3825331 | 0009910452 јо. 


0.163393469 0.008566066 0.225 
0.170164244 0.007874719 0.235 


Meriodinal plane shape and basic dimensions 
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Figure: Stepwise simulation 
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Figure: Blade Geometry isometric 3D view 
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Figure: Mesh Elements in Turbo mesh 
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Figure: Performance results shown by ANSYS with 67% efficiency 
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Figure 21. Velocity Streamlines at Blade ТЕ 








Figure: Velocity Streamlines at Blade Trailing Edge 
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Figure: Performance results shown by ANSYS with 53% efficiency 


32 


EE 
a | омен _ 
ine Fon нна 
сасе | шз _ 


ENE 
эшит | oen) С 


m Бо ПЕ 
% 


Figure: Performance results shown by ANSYS with 52% efficiency 





Figure: Velocity vectors span shown in ANSYS 


33 





Figure: Performance results shown by ANSYS with 45% efficiency 
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